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Abstract 



We expose a new symmetry for linear perturbations around a solution of non-linear Fierz-Pauli massive 
gravity plus a bare cosmological constant. The cosmological constant is chosen such that the background 
^ ' metric is flat while the Stiickelberg fields have a non-trivial profile. Around this background, at linear 

D ■ order the new symmetry reduces the propagating degrees of freedom to those of General Relativity, namely 

i-C \ the massless helicity 2 modes only. We discuss the physical consequences and possible applications of these 

findings. 

> 

O ■ 1 Introduction 

00 

' One of the most important properties of general relativity (GR) is that it propagates two helicity-2 interacting 
^j' ■ degrees of freedom. This property is ensured by the gauge symmetry associated with general coordinate 
! invariance, together with the particular structure of the Einstein-Hilbert action. The fact that long-range 
CN ' fifth forces, for example due to light scalars, are not observed in the Solar system, suggests that GR is 
the correct theory of gravity near the Earth. On the other hand, serious theoretical issues with Einstein 
gravity in the infrared, as the cosmological constant problem and the need to explain the observed current 
acceleration of the universe, motivate to consider alternative theories of gravity. Typically, these theories 
involve modifications of GR at short and/or large distances, and provide compelling explanations for the 
observed acceleration of the universe. However, in most cases they spoil the properties that characterize GR, 
such as coordinate invariance or the number of degrees of freedom. These additional degrees of freedom may 
have unacceptable pathologies (for example, they may be ghosts), or lead to predictions in disagreement 
with current observations. However, as we are going to discuss here, on some backgrounds, new gauge 
symmetries might emerge, which are able to remove the new undesired degrees of freedom. 

One of the simplest modification of GR consists of adding a mass term for the graviton to the Einstein- 
Hilbert (EH) action, leading to a theory of massive gravity. This possibility was proposed long time ago by 
Fierz and Pauli [I], who considered a mass term which is uniquely defined at linear order in perturbations 
around flat space. Here we consider a theory that minimally extends the Fierz-Pauli (FP) action beyond 
linear order in perturbations which is given by (we set M 2 /2 = 1): 

I = I EH + I FP = J d^x^g- [R - 2A] - ^ J d 4 x^h^h p(7 {g^g™ - gTsT) ■ (1) 



In the previous formula is the dynamical metric in four dimensions, and = g^ u — rj^ represents 
the displacement of the dynamical metric from a fiducial, fixed metric that we choose to be the Minkowski 
spacetime. Notice that we have also added a bare cosmological constant term to the usual Einstein-Hilbert 
piece, controlled by the parameter A. 

In the case A = 0, Fierz and Pauli proved that the action describing quadratic perturbations around 
Minkowski space contains five massless degrees of freedom: two tensor, two vectors, and one scalar pQ. 
Some time later, Boulware and Deser showed that an additional propagating degree of freedom, a ghost, 
arises when a point-like source is added to flat space, whose kinetic term is weighted by the inverse of the 
source mass [2]. In total, around a generic background this theory propagates six degrees of freedom. This 
is expected since general coordinate invariance is broken by the mass term. See [3] for a recent review. 

However, it is interesting to ask whether a background solution for the theory ([T]) exist, around which 
a smaller number of degrees of freedom propagate due to the emergence of new gauge symmetries. The 
answer is affirmative, as we will show in this work. Apart from the Minkowski spacetime, other exact 
vacuum solutions have been found for the action ([I]) with no cosmological constant [1]. In this paper, we 
first generalize these solutions to include an arbitrary bare cosmological constant A, and then analyse in 
detail the dynamics of linear perturbations around a particularly interesting and simple configuration, in 
which the effect of the graviton mass and the cosmological constant terms compensate each other, and lead 
to a flat physical background. Around this configuration, a new gauge symmetry emerges in the quadratic 
action for perturbations, which removes the dynamical vector and scalar degrees of freedom, leaving the 
massless helicity-2 mode as the only dynamical state. Consequently, the dynamics of linearized fluctuations 
for this theory around our flat background behaves exactly as in GR. 

Our findings explicitly show that enhanced symmetry points exist in the space of background solutions 
of modified gravity models, such as massive gravity. These configurations might represent examples of 
consistent backgrounds where we can live on, and in which the infrared issues of GR could be addressed 
more successfully than in GR. 

The paper is organized as follows. In section 2, we discuss self-accelerating solutions in the generalised 
Fierz-Pauli model, including a bare cosmological constant. We show that, imposing a particular tuning 
of parameters, the physical metric reduces to the Minkowski spacetime. In section 3, we study linear 
perturbations around this non-trivial flat space solution, using the ADM formalism. We split the analysis 
in tensor, vector and scalar modes with respect to the spatial metric and show that it only propagates the 
massless tensor modes as in GR. In section 4, we explicitly construct the gauge symmetry which reduces 
the number of physical degrees of freedom around our specific background, and we show how to express it 
in a covariant way. Finally, we conclude in section 5 with a discussion of our results. 

2 New exact solutions 

In order to study the theory (pQ) and its exact solutions, it is convenient to implement the Stiickelberg trick 
to recover diffeomorphism invariance following [5j. This can be achieved by introducing a covariantization 
of the metric displacement in terms of the following definition 



Therefore, the quantity H^ u is given by a combination of the physical metric g^ u , and a metric in field space 



9fiu = + 




Va/3 d^dvfj) 1 



(3) 



The FP part of the Lagrangian density, in terms of the quantity , can be written as 




(4) 
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The Stiickelberg fields 4> a have been introduced to restore reparametrisation invariance, so that the previous 
action is invariant under x a — > x a + £ a , provided these Stiickelberg fields 4> a transform as scalars with 
respect to space-time symmetries (their indexes are raised and lowered by means of the fiducial metric rj a ^). 
Choosing the unitary gauge, defined by (fr a = x a , one obtains H^ u = h^ u which leads to the theory written 
as in eq. ([T|). The simplest background solution for this theory, in absence of bare cosmological constant, is 
Minkowski space, given by 

9av = r,^ , <p = a? , H IU , = , A = 0. (5) 

Salam and Strathdee found other vacuum solutions, that we now generalize to the case of a non-zero 
cosmological constant (see Appendix [A] for a full derivation). We write them in a coordinate system that 
is particularly useful for studying cosmology, where the metric reduces to the flat slicing of a maximally 
symmetric space, namely 

ds 2 = -dt 2 + e 2M (dx\ + dx\ + dxfj , (6) 
while the Stiickelberg fields acquire the following non-trivial background profile 



[3 1 [ , /2sinhmt + rh 2 r 2 e lht \ _ ^ ,, 

\ arctanh „— ^ — =-r —rare +arctanh( 

V 2 fh [ \2 cosh fht — fh z r A e mt ) v 



mre mt ) 



(7) 



\e^x l (z = 1,2,3). (8) 



In the previous equations we dubbed for simplicity r = \J x\ + x\ + x 2 , and called 

9 m 2 A 

m 2 = — + - • (9) 

Notice that the Hubble scale is set by the FP mass and the cosmological constant, thus a solution exhibits 
the acceleration if fh > 0. In absence of cosmological constant, A = 0, and for a non-zero m, the solution 
exhibits the self-acceleration, i.e. an acceleration without cosmological constant. 

The general solution discussed above admits an interesting limit leading to flat space, when choosing the 
cosmological constant A = —(3/4) m 2 . It is easy to check that this background solution for the metric and 
Stiickelberg fields reads 

/3 1 3 

-x^ , H tiv = --r] liv , A = --m 2 . (10) 

Notice that the previous solution is different from the simplest Minkowski solution in the absence of cosmo- 
logical constant term (eqs. ([5])), by a constant in the Stiickelberg fields. The field space metric f^ v , defined 
in eq. (|3|), retains the same symmetries of the physical space-time metric g^ u . 

In the next section, we will study the dynamics of perturbations around the background of eq. ()10p . 
showing that a new emerging gauge symmetry reduces the number of degrees of freedom to the ones of pure 
GR with no mass term and no cosmological constant. 

3 Dynamics of perturbations 

3.1 Definition of perturbations 

We follow an ADM approach and use Latin indexes to characterize quantities in the three spatial dimensions. 
The metric reads 

ds 2 = —N 2 dt 2 + (dx i + N { dt) (dx j + N j dt) . (11) 
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Perturbations to the lapse function and the shift vector are defined as 

N = 1 + A , N { = B\ 
and to the three dimensional metric as 

1 



Hi 



Sij + 2/t£ 5ij + 2 ( £ij — -$ij£k I 



(12) 
(13) 



where the curvature perturbation is defined as 1Z = hi — \£^- Using the extrinsic curvature K%j 



(ViiVj + VjN{ — jij), we express the Einstein-Hilbert part of the action ([I]) in the usual form 



Ieh 



dtd 3 x Ny/j 



R-2A + K\K\ - K 2 - 2A 



(14) 



Furthermore, it is convenient to decompose the various quantities into scalar, vector and tensors with 
respect to the 3D spatial metric as 



Bi 
f- ■ 



Si + B 



■ i • 



i j - 



with Si and F{ transverse vectors, and transverse traceless tensor, namely 



Sj — F ,i — K — h\ i — 
The perturbation of Stiickelberg fields split as 

= (<P, Xi + Vv) 

with x l i = 0. 



0. 



(15) 
(16) 

(17) 

(18) 



3.2 Action quadratic in perturbations 

At linear order, the action of perturbations vanishes in agreement with (|10p being a solution of ([1]). At 
quadratic order in perturbations around the background solution (|10p . the EH contribution to the action 
(H|) reads 

jguad = J dt d 3 x | hij h ij + h ij Ahij + 2hhij h ij + ^(S i - (Si - fj) . 

-6TZ 2 + ATI A(B - E) - AA ATI + 2 d{RJ&Tl 



( 3 1 

-2A ( 3AK + -K 2 + (A + 1Z)AE + -(AE) 2 - E^E' 13 - F^F 1 ' 3 



V 2 

where dots represent time derivatives. The FP part is given by 

,2 



) } , (19) 



jquad 
l FP 



m 

"T 



dtd 3 x 



U> 3\/~6A + 3\/~6AE + 9V&R, -6Atp- 3ip 
9 



+A 



-A — 3 AE — 91Z + 3^/6 Atp 



63 

-21 AE - y7£ + 9\/6A^ 



-3 ( A^- X I-AE 



6hijh ij - GFijF 1 ' 3 ' - QE tij E' i3 + 3(AE) 2 } . (20) 
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Adding these two pieces, and setting A = —3/ 4m 2 , we obtain that the tensor and vector contributions 
organize in such a way to become exactly the ones of pure GR in the Minkowski background, namely 



quad 
tens,vect 



dtd 6 x 



hijh 1 ' + h^Ah ij + -[S l - F 



Si - R 



(21) 



Therefore, they behave as in pure GR, describing the propagation of helicity-2 states. One can prove this 
standard fact by a Hamiltonian analysis of the action f)21 [) . as shown in detail in Appendix [Bj 

After fixing A = — 3/4m 2 , the Lagrangian density associated with the total scalar contribution, up to a 
total derivative, is 



C 



quad 
seal 



{ - 67t 2 + 4TZ A(B — E) — 4A A1Z - 2KA1Z} 



m 

T 1 



3V6A + 3V6AE + 9V6K -6A?p-3ip 



-A 



-TZ 



-A - 9 AE - 27K + 3\/6 Atp 



81 

-27 AE- — 1l + 9VGAiP 



3 Aip 



AE 



}■ 



(22) 



3.3 Hamiltonian analysis of the scalar sector 

It is convenient to perform a Hamiltonian analysis of this system, in order to count the number of dynamical 
degrees of freedom in the scalar sector described by the quadratic Lagrangian ([22]) . We first notice that this 
Lagrangian does not contain time derivatives of A, B, and ip; these quantities will then be associated with 
constraints. Therefore, we have in principle three dynamical scalar degrees of freedom: TZ, E and (p. The 
conjugate momenta associated with these dynamical variables are 

gr 

U n = — £ = AAB - 4AE - Y2TZ , (23) 

dn 

U E = ^ = -4ATZ, (24) 
d£ 

IF = ^ = --m 2 (V6A + VQAE-2^-2A^ + 3V6Tl). (25) 
dip 4 \ / 

The scalar Hamiltonian is defined as 

H s = TZU n + EU £ + ipU' fi -Cs + ACa + BCb + ipC^ (26) 

where we introduce the constraint Ca = dC s /dA, and similarly for B and ip, so that Ti s is written in the 
following way 

H s = ^ 9 (A~ 1 n B ) 2 -6 (A- 1 ^) n^ + 4n^ (^- + 9VGn\ + 12 A ( + 4ft 2 ) .(27) 



The constraints read 



C A = -4A^-j|lF, (28) 

C B = -n £ , (29) 
C$ = AW. (30) 
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These are first class constraints, since their mutual Poisson brackets vanish: 



{C A , C B } = {C A , C^} = {C 4 „ C B } = 0. (31) 
Moreover, the Poisson brackets among the constraints and the Hamiltonian satisfy the following relations 

{c A , h s } = -C B , (32) 
{C^Us} = 0, (33) 

{C B ,n s } = ^1^, (34) 

which imply that the constraints are preserved under time evolution. After imposing these constraints, it 
is straightforward to check that the Hamiltonian vanishes. Indeed three first class constraints are able to 
remove all the phase-space dynamical degrees of freedom. Consequently, also the scalar sector does not 
contain any physical degree of freedom. 

To summarize, quadratic perturbations around the flat-space background (|10p lead to the propagation 
of only helicity-2 modes, exactly as in GR. 



3.4 The new gauge symmetry 

We can understand the result of the previous section in terms of gauge symmetries holding for the action of 
quadratic perturbations around the background solution (|10p . Starting from our Hamiltonian analysis, we 
can implement the standard rules for obtaining the gauge symmetry associated with the first class constraints 
(see for example [B]). The infinitesimal transformations that leave the quadratic scalar Lagrangian (|22|) 
invariant, up to total derivatives, read 



dip = 




(35) 


8ip = 


h + yf e B , 


(36) 


STZ = 


o, 


(37) 


5E = 




(38) 


5A = 


£A , 


(39) 


SB = 


£B - £A , 


(40) 



for three arbitrary functions ca, £b and e^. The parameters €a and es are associated with the diffeomorphism 
invariance of the FP action, once the Stiickelberg fields are introduced. The symmetry associated with ew,, 
that acts only on the Stiickelberg scalars, is the new symmetry. 

It is not difficult to express these symmetries in a covariant way. Adding the cosmological constant 
contribution to the covariant FP action tuning A = — 3/4 m 2 , and expanding at quadratic order in 
perturbations around our solution (|10p . one finds the following covariant Lagrangian density 

A =^ 2 (^-\/|^) 2 - W 
It is easy to check that this Lagrangian is invariant, at quadratic order, under the transformation 

5^ = 5P - + X 11 with = 0, (42) 

hfj, v = h^y - . (43) 
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Focussing on scalar perturbations, the symmetry (|42p corresponds exactly to the extra symmetry parame- 
terized by the function in eqs ([55]) and The existence of this gauge symmetry, besides the standard 
diffeomorphism invariance, implies that the scalar and vector degrees of freedom do not propagate around 
the background (|l(jp . Note that this resulting quadratic action (|4ip is exactly the same as a graviton mass 
term proportional to (H^) 2 . It is known that this mass term does not modify GR, i.e. it leads only to a 
massless transverse-traceless graviton jjj, which is consistent with our result. On a different setting, it was 
also shown in [8] that, thanks to a new gauge symmetry, massive gravity on a de Sitter background does 
not propagate scalar zero modes, when the size of the cosmological constant is tuned to a precise value 
depending on the graviton mass. 

4 Discussion 

We constructed new solutions for a covariantized non-linear Fierz-Pauli theory of massive gravity, equipped 
with a cosmological constant A. We focused on a particularly simple configuration corresponding to 
Minkowski space with a non-trivial profile for the Stiickelberg fields. The space-time flatness is achieved by 
tuning A to a particular value related to the graviton mass. We then studied linear perturbations around 
this configuration, and showed that it only propagates the massless helicity-2 modes, in the same way as 
pure GR without a cosmological constant. We interpreted this behavior as due to a new gauge symmetry 
acting on the Stiickelberg fields, in addition to diffeomorphism invariance. 

It would be important to understand whether this new symmetry is only associated with the particular 
set-up discussed here, or whether it arises more generally. In other non-linear massive gravity models, 
several examples of self-accelerating solutions have been found [UJ [TUJ, [TTJ [JSJ EH EH IS]- Recently, [H] 
found similar situations where there are no propagating scalar and vector modes around self-accelerating 
open-FRW solutions in the non-linear massive gravity models proposed by de Rham, Gabadadze and Tolley 
|16j . See also [15] for examples of self-accelerating solutions in non-linear theories of massive gravity, around 
which there are no additional degrees of freedom besides tensor modes. One one hand, a common feature 
between our solution and the set-up of [2] is that the metric in field space preserves the same symmetries 
as the physical metric, which in our case correspond to the fact that both are Minkowski spacetime. On 
the other hand, there are other self-accelerating de Sitter solutions in the same theory for which the field 
space metric is not invariant under same isometries as the physical metric [9j \10\ [TT] . Those solutions are 
indeed similar to the ones obtained in this article, but where the cosmological constant is not tuned with 
the mass term to give flat space. The non-trivial profiles for the Stiickelberg fields render the analysis of the 
dynamics of fluctuations more difficult. However, it would be interesting to study whether new symmetries 
also arise while considering perturbations around these backgrounds. 

Finally, it is also important to explore whether the new gauge symmetry holds only at linear order in 
perturbations, or survive at higher orders. If the latter is the case, then these configurations may represent 
examples of consistent backgrounds where we can live on, and in which the infrared issues of GR could be 
addressed more successfully than in GR. For example, one could envisage a mechanism that dynamically 
tunes the graviton mass with the cosmological constant, leading to flat space also in the presence of A. We 
will investigate these interesting questions elsewhere. 
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A Derivation of the exact solutions 



rj^. Since we are interested in static 

(44) 



We begin by choosing the gauge <^ = x^, which implies that f^ v 
spherical solutions to the action ([1]), we introduce the ansatz 

ds 2 = -C(r) dt 2 + A(r) dr 2 + 2D(r) dtdr + B(r)dQ 2 , 

where dU 2 = d6 2 + sivPOdcf) 2 , and to simplify expressions, one can further choose to write the field space 
metric fav in spherical coordinates as f^dx^dx" = —dt 2 + dr 2 + r 2 d£l 2 . 

After introducing the ansatz (|44p in the equations of motion derived from the action (pQ), one can show 
there is a constraint which leads to two branches of solutions: one with D(r) = and the other where 
B{r) = 3r 2 /2. The first branch was studied in |17j . but here, we are interested in the second branch, where 
exact solutions were initially found by Salam and Strathdee [I]. In this second class, a new constraint 
enforces A(r) = A(r)C(r) + D 2 (r) = Aq = const, and the rest of equations of motion admit the following 
solution 



A(r) 
C(r) 



3A 



where 



2 
3A 



p(r) 



(p(r) + a + l), 
(1 -p(r)), 



B(r) 



—r 
2 



(45) 



D ( r ) = 2 A oVp( r )(p( r ) + 



c 
r 



2m 2 r 2 



3A 



3/2 ' 



a 



9 Ac 



(46) 



and rh is defined as in (|9|). Ao and c are two integration constants that obey c > and < \^Aq < 2/3 for 
D(r) to be real. In this work we only consider c = and the extremal value Ao = 4/9. For these values, the 
solution (|45p can be re-casted in the simple FRW form of eq. ([6]) by the following coordinate transformation 



3 1_ 

2 rh 

| e mt x i 



arctanh 



( 



2 sinh fht + rh 2 r 2 e mi 



\ 2 cosh fht 
(» = 1,2,3). 



mr e 



■nit 



arctanh (rh r e mt ) 



(47) 



The last transformation implies that r 



3 „fht 



e r, while 9 and <f> remain unchanged. Furthermore, the 



transformations change our initial gauge so that / MJ/ ^ rj^, since the initial Stiickelberg fields ^ = x M get 
modified to those in (J7J). 



B Hamiltonian analysis of tensor and vector degrees of freedom 



The canonical momenta associated to the tensor and vector modes in the action (12111 are 



IP, 



5C 



quad 
tens,vect 



5C 



Shij 

quad 
tens,vect 



2h i] 



A(S l 



(48) 



where ^tensvect ^ associated Lagrangian density of the action ([2T]) . In momentum's language, the Lagrangian 



C tens,vect then reads 



jquad 

tens,vect 



dtd 6 x 



iqkj + YVpFt -Ht-Hv- S t U l 



(49) 
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where 

■Ht = injfnj - fjiAhij, n v = -inf A" 1 ^, (50) 

with A -1 the inverse Laplace operator. Let us first consider the tensor modes, which do not have any 
associated constraint. Therefore, both tensor modes are physical degrees of freedom, and their Hamiltonian 
is exactly that of GR (see for example [IB]). In the vector case, F l is the only dynamical degree, while S l 
appears as a Lagrange multiplier, enforcing the constraints IT^ = 0. These constraints are first class since 
they commute among each other and with the Hamiltonian T~L V , thus they represent the two gauge modes 
associated with vector perturbations. In summary, the tensor modes as the same as in GR and there are no 
physical vector degrees of freedom. 
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